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SUMMARY 

The  objective  of  this  program  on  coherent  optical  adaptive  techniques 
(COAT)  is  to  investigate  algorithms  and  techniques  that  can  reduce  the  beam 
distortions  caused  by  thermal  blooming.  This  report  covers  the  contract 
period  from  1 November  1975  through  31  Mary  1976.  The  work  during  this 
period  has  centered  on  the  construction  and  preliminary  testing  of  a deform- 
able mirror  and  associated  electronics. 

The  construction  of  a 37-element,  all-beryllium  deformable  mirror 
has  been  completed.  Final  polishing  of  the  mirror  is  still  not  complete 
because  of  the  inability  of  a supplier  to  deliver  the  bulk  of  the  37  piezoelectric 
stack  driver  elements.  Only  1 9 of  the  37  stacks  have  been  installed  at  this 
writing.  Some  preliminary  tests  on  the  mirror  hav~  determined  its  excur- 
sion sensitivity,  influence  function,  and  frequency  response.  At  present, 
the  mirror  delivers  only  0.28  |jim  of  surface  motion  for  1 50  V of  drive, 
slightly  less  than  one  half  the  design  goal.  The  frequency  response  shows 

some  unusual  resonance  behavior,  but  is  usable  to  at  least  30  kHz.  The 

1 . 5 

influence  function  has  a shape  given  by  IN(r)  = exp(-Br  ’ ) and  has  a value 
of  IN  = 0.23  at  the  center  of  the  adjacent  actuator.  The  faceplate  will  be 
cut  thinner  to  reduce  this  value  to  about  0.  15. 

The  computer  simulation  studies  of  irradiance  tailoring  reported  in 
the  previous  report  were  continued  to  correct  an  error  in  the  results  pre- 
sented there.  The  basic  conclusion  is  unchanged,  however:  the  best  way  to 
reduce  thermal  blooming  is  to  increase  the  transmit  aperture  size  and 
make  the  beam  as  uniform  as  possible. 

We  have  completed  construction  of  high-voltage  mirror  drivers  and 
a dither  injector  for  use  with  the  deformable  mirror.  Construction  has 
begun  on  a 14-function  Zernike  polynomial  generator.  Computer  simulation 
studies  have  identified  which  Zernike  polynomials  make  the  greatest  con- 
tribution to  blooming  compensation.  A computer  code  has  been  developed 
that  provides  data  on  generating  Zernike-polynomial  mode  shapes  with  a 
finite  number  of  deformable  mirror  elements  with  minimum  mean  square 
error  between  the  actual  Zernike  shape  and  the  mirror-produced  surface. 


During  the  next  period  the  deformable  mirror  construction  and 
characterization  will  be  completed.  The  Zernike-polynomial  generator 
will  be  completed,  and  performance  tests  will  begin  with  the  deformable 
mirror  in  the  DARPA/RADC  COAT  system. 


8 


I.  INTRODUCTION 

A.  Program  Objectives 

The  primary  objective  of  this  program  on  coherent  optical  adaptive 
techniques  (COAT)  is  to  analyze  and  experimentally  demonstrate  adaptive 
multidither  correction  algorithms  that  can  reduce  beam  distortions  caused 
by  thermal  blooming.  The  use  of  fixed  transmitter  intensity  profiles  for 
reducing  thermal  blooming  will  also  be  investigated. 

B.  Research  Program  Plan 

The  research  program  utilizes  the  21 -channel  DARPA/RADC  experi- 
mental COAT  system  built  and  tested  on  contracts  F30602-73-C-0248  and 
F30602-75-C  -0001 . Computer  simulation  codes  developed  on  these  con- 
tracts and  on  other  programs  (e.  g.  , NSWC  contract  N60921 -74-C-0249) 
have  been  used  for  the  analytical  portions  of  this  contract.  These  ''odes 
model  the  operation  of  several  types  of  COAT  servomechanisms  as  well  as 
the  time-dependent  propagation  of  optical  beams  in  an  absorbing  and  turbu- 
lent medium.  The  experimental  investigations  require  construction  of  a 
deformable  mirror  for  the  COAT  system  as  part  of  the  program. 

To  accomplish  the  objectives  of  this  contract,  a 12-month  research 
program  consisting  of  three  major  tasks  was  developed.  Task  1 is  an 
analytical  task  and  is  largely  completed.  Task  2 provides  for  the  construc- 
tion of  a deformable  mirror  and  is  also  nearly  complete.  In  Task  3 we  will 
use  this  mirror  to  study  thermal  blooming  compensation  with  zonal  multi- 
dither control  and  with  modal  control  using  Zernike-polynomial  modes. 

Amendment  No.  1 to  the  contract  has  provided  funds  to  characterize 
the  beryllium  deformable  mirror  being  built  and  to  assess  the  applicability 
of  this  design  to  high-power  laser  COAT  systems.  The  characterization  of 
the  mirror  will  occur  prior  to  the  start  of  Task  3. 


II. 


TECHNICAL  ACCOMPLISHMENTS 


A.  Zernike  Polynomial  COAT  (ZEP-COAT) 

At  the  beginning  of  this  program,  we  had  some  hope  that  multidither 
control  of  Zernike  polynomial  functions  (modal  control)  might  produce  a 
greater  degree  of  compensation  for  thermal  blooming  than  localized  (zonal) 
control  of  the  transmitted  phasefront.  As  a result  of  our  computer  simula- 
tion work  and  experimental  results  at  Hughes  and  at  Lincoln  Laboratory, 
and  through  a better  understanding  of  the  limitations  of  phase  compensation 
of  thermal  blooming,  we  no  longer  believe  that  modal  control  has  any 
inherent  advantage  over  zonal  control  in  the  degree  of  compensation  that  it 
can  produce.  This  conclusion  is  certainly  true  for  turbulence  compensa- 
tion using  Zernike  polynomial  modal  control,  as  shown  by  the  results  in 
Figure  1.  For  turbulence  compensation,  it  is  the  number  of  control  chan- 
nels (number  of  degrees  of  freedom)  that  counts,  not  what  functions  are 
controlled. 

For  thermal  blooming,  however,  it  may  be  possible  to  obtain  sig- 
nificant blooming  correction  with  far  fewer  control  channels  than  with  a 
zonal  control  system.  The  results  of  Figure  1 also  indicate  that  signifi- 
cant turbulence  correction  is  possible  with  a limited  number  of  Zernike 
control  channels.  With  fewer  channels  the  complexity  of  the  control  sys- 
tem is  reduced,  and  faster  servo  response  is  possible  when  the  dither 
mirror  response  limits  the  servo  bandwidth.  This  last  conclusion  results 
from  the  computer-determined  result  that  the  maximum  servo  bandwidth, 
(fs)max»  f°r  a number  of  sine/cosine  dither  channels,  Nc>  is  given 

by 
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where  f is  the  maximum  dither  frequency  allowed  by  the  deformable 
max 

dither  mirror.  Note  that  N in  Eq.  (1)  is  the  number  of  control  channels, 
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Analytical  comparison  of  turbulence  compensation  with 
Zernike  polynomials  and  zonal  COAT  control.  The 
transmitter  diameter  is  D_ . 


not  the  number  of  dither  frequencies  (one  half  of  Nc  in  a sine/cosine 
system),  and  N is  assumed  to  be  a multiple  of  2. 

In  principle,  a multidither  COAT  system  can  operate  with  any 
function  set  that  is  orthogonal  (or  nearly  so)  over  the  transmitting  aper- 
ture. That  is,  the  correction  phase,  <t>c»  can  be  generated  according  to 
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where  the  two-dimensional  functions  f are  orthogonal  over  the  correction 

aperture,  a^  are  the  correction  amplitudes,  and  and  are  the  dither 

amplitude  and  frequency  for  the  n*b  servo  channel.  For  a zonal  system, 

the  functions  f are  the  influence  functions  for  the  N deformable  mirror 
n 

actuators. 

For  optical  systems  with  circular  apertures,  a logical  choice  of  an 
alternate  orthogonal  correction  set  is  the  Zernike  polynomials.  These  func- 
tions are  orthogonal  over  a circular  aperture,  and  the  lowest-order  func- 
tions are  the  classical  aberrations  of  tilt,  focus,  coma,  etc.,  most  of  which 
are  known  to  have  large  values  in  phase  distortions  produced  by  turbulence 
and  thermal  blooming.  We  have  given  the  name  "ZEP-COAT"  to  a multi- 
dither modal  COAT  system  that  uses  Zernike  modes. 

Using  computer  simulation  codes,  we  have  determined  that  only  the 
7 polynomials  indicated  with  an  asterisk  in  Table  I are  necessary  for  good 
correction  in  moderate  thermal  blooming.  In  determining  that  these 
7 polynomials  were  the  significant  ones,  a moderate  blooming  case  was 

chosen  (P  . ~ P ,).  The  desired  cw  correction  phasefront  was 

' transmit  critical  a ^ 

calculated  using  a function-maximization  routine.  * This  phasefront  was 


PW.  P.  Brown,  Jr.,  "Computer  Simulation  of  Adaptive  Optical  Systems,  " 
Final  Report  on  Contract  N60921  -74-C-0249,  September,  1 975. 

bJ.  E.  Pearson,  et  al. , "High  Power,  Closed  Loop,  Adaptive  System 
(HICLAS)  Study,"  Final  Report  on  Contract  N60921  -76-C-0008. 


Table  1.  Zernike  Polynomials  Used  in  ZEP-COAT  Simulation 


n 

Z 

n 

Aberration  Name 

** 

Pn 

*1 

X 

Horizontal  tilt 

2 

2a 

Y 

Vertical  tilt 

2 

*3 

2R2  - 1 

Refocus 

1.5 

*4 

2 2 
X c - Yc 

Astigmatism 

3 

5a 

XY 

Astigmatism 

3 

* 6 

X(3R2  - 2) 

Coma  (Y-symmetry) 

4 

7a 

Y(3R2  - 2) 

Coma  (X-symmetry) 

4 

8 

X(X2  - 3Y2) 

120  degree  (5th  order, 
Y - symmetry) 

4 

9a 

Y(Y2  - 3X2) 

120  degree  (5th  order, 
X-symmetry) 

4 

*10 

6R4  - 6R2  + 1 

Spherical  aberration 

2.  5 

1 1 

(X2  - Y2)(4R2  - 3) 

Astigmatism  (5th  order) 

5 

12 

XY  (4R2  - 3) 

Astigmatism  (5th  order) 

5 

* 

13 

X(10R4  - 12R2  + 3) 

5th  order  coma  (Y-symmetry) 

6 

14a 

Y ( 1 OR4  - 12K2  + 3) 

5th  order  coma  (X-symmetry) 

6 

*1  5 

20R6  - 3 OR4  + 12R2  - 1 

5th-order  spherical  aberration 

3.  5 

where 

X = — 7=  = R cos  0,  Y = y = R sin  0,  and  R = x2  + y2  s 1 
/2  /2 

cl 

These  five  polynomials  are  not  used  for  compensating  only  thermal 

blooming  because  of  the  symmetries  in  the  blooming  distortions;  the 
transverse  wind  is  in  the  X-direction  (see  Ref.  11).  They  are  used  for 

turbulence  compensation. 

^Principal  polynomials  needed  for  blooming  compensation. 

jJcjJc 

Relative  dither  amplitude  for  constant  effective  modulation  index. 

T1994 
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then  fitted  by  mean3  of  an  expansion  using  the  first  44  Zernike  polynomials. 
Table  2 lists  the  results  of  this  calculation  and  Figure  2 shows  the  uncor- 
rected and  corrected  beam  profiles  calculated  for  the  conditions  shown  in 
Table  3.  Although  there  are  some  contributions  from  higher  order  poly- 
nomials, their  amplitudes  are  either  le3s  than  1/10  of  the  smallest  ampli- 
tude among  the  7 polynomials  noted  in  Tables  1 and  2,  or  are  of  too  high  an 
order  to  be  produced  by  a 37-element  mirror.  For  the  planned  experi- 
ments, however,  we  will  implement  all  14  polynomials.  The  ZEP-COAT 
system  can  thus  be  used  for  turbulence  compensation  studies  if  so  desired. 

Implementing  a ZEP-COAT  capability  into  the  DARPA/RADC  hard- 
ware is  a straightforward  procedure,  as  indicated  in  Figure  3.  A "poly- 
nomial generator"  is  merely  inserted  between  the  control  channels 
(synchronous  detectors  plus  low-pass  filters)  and  the  drivers  for  the 
deformable  mirror.  The  polynomial  generator  is  composed  of  several 
groups  of  weighting  matrices  that  distribute  the  error  signal  across  the 
mirror  in  the  correct  manner  to  produce  the  desired  polynomial  spatial 
mode.  Since  the  dither  and  correction  signals  will  be  put  onto  a single 
mirror  for  our  experiments,  the  low-voltage  dither  and  correction  signals 
will  be  combined,  put  through  the  polynomial  separated  and  separately 

amplified,  and  then  recombined  onto  the  mirror.  The  process  is  illustrated 
in  Figure  4. 

A 14-channel  polynomial  generator  is  nearly  complete;  37  mirror 
actuator  drivers,  each  capable  of  ±150  V output,  and  a dither/correction 
combiner  (dither  injector)  have  already  been  completed.  The  interconnec- 
tion of  these  units  for  either  a zonal  or  a modal  control  system  is  shown  in 
Figure  5.  Because  of  the  excessive  delays  in  the  delivery  of  the  piezo- 
electric transducer  (PZT)  material  for  the  deformable  mirror  (Section  B 
of  this  chapter),  the  technical  portion  of  the  contract  will  be  extended  in 
time,  but  at  no  increase  in  cost,  until  the  middle  of  November  1976  in 
order  to  complete  the  necessary  thermal  blooming  studies. 
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Figure  2.  Thermal  blooming  computer  simulation  data  for 
conditions  listed  in  Table  3.  (a)  Uncorrected 

focal-plane  beam  profile:  Strehl  ratio  = 0.  11. 

(b)  Corrected  beam  using  function-maximization 
routine:  Strehl  raUo  = 0.  34.  The  phase  correc- 
tion used  here  is  used  to  generate  the  Zernike 
polynomial  amplitudes  in  Table  2. 


Table  2.  Expansion  of  Thermal  Blooming  Phase  Correction  in  Terms  of 
Zernike  Polynomials  Rm(p)  • sin  m6  on  Rm(p)  • cos  mfl 


Numbe 


Spherical  aberrati 


oma  (5th  order! 


To  be  implemented 


Table  3.  Assumed  Conditions  for  Thermal  Blooming  Calculations 


— x 

Range:  2 km  Absorption  coefficient:  4x10  m 

Wavelength:  10.6  pm  Slew  rate:  10  mrad/sec 

Transmitter  diameter:  1 m,  Transverse  wind  velocity:  2 m/sec  in  +X 

uniform  beam  direction 

Resultant  uncorrected  Strehl  ratio:  0.  39 
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The  computer  simulations  being  used  to  analyze  ZEP-COAT 
performance  were  discussed  in  the  previous  contract  report.  We  also 
reported  that  we  were  observing  an  instability  in  the  control  system.  We 
have  determined  that  the  instability  was  caused  by  an  improper  setting  of 
the  gain  and  modulation  index  for  the  channels.  We  have  now  optimized  both 
these  quantities  and  the  ZEP-COAT  servo  system  exhibits  stable  behavior 
under  all  conditions.  If  the  dither  signals  applied  to  the  mirror  are  repre- 
sented by 


^(pn  sinWnt)  Zn(r,  6) 
n 


(3) 


where  Z (r,  6)  are  the  Zernike  polynomials,  the  the  relative  values  of  (3 
n n 

for  a 5%  effective  modulation  index  are  given  in  the  last  column  of  Table  1. 
The  values  of  Pn  depend  on  the  order  of  the  polynomial  according  to 


Pn 


<x 


n + 1 
2 


n + 1 


(4) 


where  m and  n are  given  in  Table  2. 
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The  computer  simulations  have  also  been  used  to  investigate 
multiple  convergence  states  for  ZEP-COAT  systems,  that  is,  the  modal 
equivalents  of  2Nir  ambiguities  in  zonal  systems. 

The  2Nir  ambiguity  problem  is  normally  associated  only  with  zonal 
multidither  systems.  An  analogous  problem  can  occur  with  a modal  ZEP- 
COAT  system,  however.  That  is,  there  exist  local  maxima  states  where 
the  dither  control  system  can  stop  (obtain  a zero  error  signal),  but  which 
are  nonoptimum  for  maximum  focal-plane  ir radiance.  Taken  individually, 
the  two  most  obvious  cases  are  tilt  (the  servo  can  place  a sidelobe  of  the 
Airy  pattern  on  a glint)  and  focus. 

The  behavior  of  the  focus  control  can  be  understood  by  referring 
to  Figure  6.  Figure  6 illustrates  the  w.'dl-known  behavior  of  the  on-axis 
intensity  of  a focused  beam:  as  the  focal  point  is  approached  from  either 
direction,  the  on-axis  intensity  goes  through  several  maxima  and  minima. 
Each  maxima  away  from  the  focal  plane  corresponds  to  a different  Fresnel 
zone,  or  a 2Nir  phase  difference  between  rays  originating  from  the  center 
of  the  focusing  lens  and  the  outer  edge.  If  the  initial  defocus  is  at  a point 
such  as  A in  Figure  6,  the  ZEP-COAT  focus  control  will  drive  the  focus 
toward  the  maximum  on-axis  intensity.  If  the  starting  point  is  at  B,  how- 
ever, the  multidither  servo  moves  toward  the  nearest  maximum  at  C and 
a nonoptimum  target  intensity  results. 

The  problem  can  also  exist  for  other  Zernike  functions.  Table  4 
lists  the  first  13  polynomials  along  with  the  location  of  the  first  secondary 
maximum  and  the  on-axis  irradiance  at  that  maximum.  At  this  point,  we 
have  only  pointed  out  the  potential  for  a servo-ambiguity  for  ZEP-COAT. 
Whether  or  not  such  ambiguities  will  occur  in  actual  systems  when  they 
attempt  to  correct  for  arbitrary  or  random  errors  is  a subject  for  future 
study. 

We  are  currently  setting  up  a computer  code  which  will  determine 
the  relative  actuator  drives  necessary  to  produce  a given  Zernike  poly- 
nomial with  a deformable  mirror.  This  code  will  be  used  to  set  the  values 
of  the  various  resistors  in  the  weighting  matrices  of  the  Zernike  polynomial 


Table  4.  2Nir  Problem  with  Zernike  Polynomials 


Zernike  Polynomial  (Aberration) 

Approximate 
Location  of 
Secondary 
Maximum, 
p.  (radians) 

On-Axis  Intensity 
at  Secondary 
Maximum  Relative 
to  that  of 

Unaberrated  Focus,  I 
' o 

1.  v|i  = P^(2p2  - 1)  [defocus] 

mam 

0.05  I 

o 

2 . 4>  = P2(6p4  - 6p2  + 1) 

5.75 

0.12  I 

o 

[spherical  aberration] 

3.  4>  = p3(20p3  - 3p4  + 12p2  - 1) 

8.  7 

0.  il  I 

o 

4.  4>  = P4X  5.  4)  = PgY 

5.  25 

0.017  I 

o 

[tilt] 

6.  4>  = P6X(3p2  - 2) 

Tiny  one  at 
7.  5 

0.00012  I 

o 

[coma] 

7.  4>  = P6Y(3 p2  - 2) 

Next  at  12 

0. 022  I 

o 

8.  4>  = P8X(10p4  - 1 2p2  + 3) 

9.  4>  = P9Y(10p4  - 1 2p2  + 3) 

14 

0.012  I 

o 

10.  4>  = P1Q(X2  - Y2) 

8.  5 

0.022  I 

o 

[astigmatism] 

11.  * = Pn(X2  - Y2)(4p2  - 3) 

13.  5 

0.03  I 

o 

12.  4>  = P12X(X2  - 3Y2) 

13.  4>  = P13Y(Y2  - 3X2) 

None  for 
p.  < 25  rad 

— 

T1997 
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generator.  The  best  fit  to  a particular  Z^  is  a function  of  the  mirror 
influence  function,  the  number  of  actuators,  and  the  interactuator  spacing. 
Proper  representation  of  a is  necessary  in  order  to  maintain  the 
orthogonality  of  the  COAT-controlled  functions.  If  the  mirror  cannot  accu- 
rately produce  the  Zernike  functions,  the  functions  that  are  produced  will 
not  be  orthogonal  and  the  COAT  control  channels  will  be  coupled  together. 
The  amount  of  the  coupling  will  be  a function  of  the  rms  difference  between 
the  mirror-produced  function  and  the  desired  Z^. 

Another  potential  source  of  coupling  with  a modal  system  such  as 
ZEP-COAT  occurs  through  the  amplitude  distribution  across  the  laser 
transmitting  aperture.  The  polynomials  listed  in  Table  1 are  orthogonal 
over  a uniform  unit  circle.  For  a nonuniform  intensity  distribution,  they 
will  not  be  orthogonal.  If  the  intensity  distribution  is  known,  a new  poly- 
nomial set,  related  to  the  Zernike  set,  can  be  found  that  is  orthogonal.  If 
the  polynomials  are  not  orthogonal,  the  result  is  again  a coupling  of  the 
COAT  servo  channels. 

Just  how  much  channel-channel  coupling  a ZEP-COAT  system  can 

2 3 

tolerate  is  not  known.  Previous  studies  of  zonal  systems  ’ have  shown 
that  20  to  30  percent  coupling  between  channels  can  be  tolerated  in  a sine/ 


cosine  dither  system  without  noticeable  degradation  of  performance.  To 
our  knowledge,  this  important  problem  of  channel-channel  coupling  in  a 
modal  control  system  has  not  been  satisfactorily  investigated. 

B.  Deformable  Mirror 

For  the  experimental  studies  with  ZEP-COAT,  a 37-actuator 
deformable  mirror  is  necessary,  and  one  is  being  built  to  replace  the 
beam-splitter/phase-shifter  array  (called  the  "phasor  matrix")  that  has 
been  successfully  used  up  to  now  with  the  DARPA/RADC  COAT  system. 

The  mirror  will  perform  the  two  functions  of  phase  dithering  (tagging)  and 
phase  correction  required  in  any  type  of  multidither  COAT  system.  Per- 
formance of  both  these  functions  with  a single  element  places  heavy  and 
conflicting  demands  on  the  devices.  The  phase  correction  function  requires 
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a large-amplitude  phase  excursion  (about  ±2  optical  wavelengths)  at 
relatively  low  frequencies  (up  to  about  1 kHz).  On  the  other  hand,  the 
dither  function  requires  low-amplitude  excursions  (±X./10)  at  much  higher 
frequencies  (10  to  30  kHz).  To  date,  no  continuous-surface,  deformable 
mirror  device  has  been  constructed  that  can  accomplish  both  of  these  func- 
tions simultaneously,  particularly  with  the  low  drive  voltages  (±125  V) 
available  from  the  DARPA/RADC  COAT  system. 

The  mirror  being  built  on  this  program  is  the  first  of  a kind,  and 
as  such  will  represent  the  state  of  the  art  for  uncooled  deformable  mirrors 
of  its  type.  Because  such  a mirror  had  not  previously  been  built,  or  even 
designed,  an  extensive  design  effort  was  initiated.  This  design  has  been 
funded  only  in  part  by  this  contract,  the  remainder  of  the  work  being  per- 
formed as  part  of  the  Hughes  Aircraft  IR&D  program. 

The  design  and  analysis  effort  started  with  Hughes  proprietary 
designs  for  multielement  deformable  mirrors  and  for  external- spring,  self- 
contained,  piezoelectrically-driven  actuator  cells.  The  basic  design  and 
pertormance  goals  are  discussed  in  an  earlier  report.* 

A schematic  of  the  mirror  is  shown  in  Figure  7.  The  entire  mirror 
body  is  made  of  beryllium.  The  37  actuators  are  placed  on  a circular  array 
with  a minimum  actuator  spacing  of  0.  5 50  in.  The  active  mirror  surface 
has  a diameter  of  3.  980  in.  The  faceplate  is  initially  0.  1 50-in. -thick 
beryllium  optical  stock,  and  is  butt-brazed  to  each  of  the  37  actuators  and 
to  the  rim  of  the  primary  backup.  The  faceplate  will  be  polished  and  coated 
with  a protective  silver  overcoat.  The  actuator  cells  are  brazed  to  the  pri- 
mary backup  at  the  same  time  as  the  faceplate  braze. 

Each  piezoelectric  (PZT)  transducer  stack  is  electrically  grounded 
at  the  mirror  faceplate  and  the  driver  voltage  is  applied  through  the  elec- 
trode at  the  other  end  of  the  PZT  stack.  With  this  arrangement,  the  trans- 
ducer stack  is  wired  internally,  thus  eliminating  the  complication  of  feedback 
wires  out  from  the  stack.  The  only  wire  to  each  actuator  is  soldered  con- 
veniently to  the  electrode  in  a manner  which  directs  the  wire  straight  out  of 
the  cell,  with  no  bends  that  require  additional  volume.  This  electrical 
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arrangement  is  compatible  with  the  electronics  of  the  DARPA/RADC 
system,  which  has  a common  ground  for  all  channel  outputs.  All  of  the 
electrical  outputs  are  brought  to  a 41 -pin  electrical  connector  that  is 
attached  to  the  dust  cover  at  the  back  of  the  mirror  (see  Figure  7). 

The  mirror  has  been  fabricated  and  is  currently  being  polished, 
although  only  18  of  the  required  37  PZT  stacks  are  presently  available.  A 
photograph  of  the  unit  is  shown  in  Figure  8.  An  analysis  performed  on  this 
mirror  predicts  that  the  lowest  structural  resonance  will  occur  at  12  to 
15  kHz,  and  appears  to  be  a shear  mode  of  the  back-plate  structure.  The 
measured  response  characteristics  for  the  central  actuator  (before  polish- 
ing) is  shown  in  Figure  9,  with  all  other  actuators  undriven,  but  with  pre- 
loads on  them  (dummy  blocks  were  loaded  into  the  outer  ring).  The 
response  curves  were  obtained  by  placing  the  mirror  in  a Twyman-Green 
interferometer  and  looking  at  the  surface  motion  of  the  central  actuator. 

The  response  of  the  other  actuators  is  similar  except  for  some  of  the 
details  of  the  fluctuations  beyond  20  kHz.  The  phase  of  the  surface  response 
is  measured  relative  to  the  sinusoidal  drive  voltage. 

The  mirror  response  is  reasonably  flat  up  to  about  20  kHz,  where  a 
rather  complicated  resonance  behavior  beings.  Using  a holographic 
interferometer,  we  conclude  that  the  fluctuations  above  20  kHz  are  a result 
of  excitation  of  complicated  actuator/faceplate  vibrational  modes.  The 
first  two  major  resonances  of  this  type  were  calculated  to  occur  at  23  kHz 
and  40.  5 kHz,  in  good  agreement  with  the  first  two  major  peaks  in  the 
amplitude  plot  in  Figure  9* 

Note  that  the  first  resonance  activity  occurs  in  the  neighborhood  of 
15  kHz,  which  is  in  good  agreement  with  the  structural  analysis.  The 
relatively  small  phase  fluctuations  indicate  that  the  mirror  is  useable  out 
to  at  least  30  kHz.  Even  a higher  frequency  response  may  be  possible  by 
filling  the  interior  of  the  unit  with  oil,  and  by  using  lower-amplitude  drive 
voltages  at  the  higher  frequencies  (as  will  be  done  for  dithers)  to  reduce 
the  drive  to  backplate  modes.  At  this  time  we  do  not  fully  understand  the 
reasons  for  the  enhanced  amplitude  response  above  20  kHz. 
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Figure  8.  Photograph  of  37-element  DARPA/RADC  deform- 
able mirror  prior  to  final  polishing. 


Figure  9.  Measured  amplitude  and  phase  response  for  central  actuator 
of  DARPA/RADC  beryllium  deformable  mirror.  Drive  to 
actuator  is  42.4  V peak,  producing  0.079  pm  of  faceplate 
motion. 


In  addition  to  the  frequency  response  measurements  shown  in 
Figure  9,  we  have  measured  the  influence  function  of  an  actuator  in  each 
of  the  3 rings  plus  the  central  actuator.  Profiles  through  the  center  of 
these  4 actuators  are  shown  in  Figure  10.  Each  actuator  can  produce 
0.28  pm  of  surface  motion  for  150  V of  drive.  As  can  be  seen,  the  influ- 
ence function  depends  on  the  distance  of  the  actuator  from  the  clamped 
edge  of  the  mirror;  the  closer  to  the  edge,  the  more  skewed  the  surface 
profile. 

The  mirror  surface  is  measured  using  a simple,  but  highly  sensi- 
tive and  accurate  technique.  A sinusoidal  drive  at  500  to  1000  Hz  is  first 
applied  to  the  appropriate  actuator  (or  group  of  actuators).  About  30  V rms 
is  used,  which  produces  a peak  surface  deflection  of  about  0.079  pm.  A 
phonograph  needle  is  then  drawn  across  the  surface  and  the  resulting  out- 
put is  observed  on  an  rms  voltmeter  as  a function  of  the  needle  location. 

The  observed  profiles  are  repeatable  to  within  a few  percent.  The  tech- 
nique is  particularly  useful  because  it  gives  the  actuator  influence  function 
independent  of  the  undriven  mirror  surface  shape.  For  use  on  polished 
mirrors,  the  needle  tip  is  dipped  in  epoxy  to  form  a nonabrasive  rounded 
surface.  Tne  epoxy  produces  very  little  loss  in  sensitivity. 

As  indicate!  in  Figure  10(a),  the  central  actuator  influence  function 
is  closely  approximated  by  the  "subgaussian"  function, 

IN(  r)  = exp(-3.  622r1- 5)  , (5) 


where  r is  the  radial  distance  in  inches  from  the  actuator  center.  Fig- 
ure 11  shows  a 3-dimensional  plot  of  this  function  and  Figure  12  is  a 
contour-line  plot  of  this  profile.  Notice  that  the  profile  is  reasonably 
symmetric,  but  is  slightly  higher  at  the  location  of  the  actuator  marked 
(a)  than  at  other  actuator  locations. 

If  we  define  a mechanical  coupling,  C , as 


C = IN(S) 
m 
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Figure  10.  Influence  function  profiles  of  beryllium  mirror. 

(a)  Center  actuator,  (b)  First  ring  actuator. 

(c)  Second  ring  actuator,  (d)  Third  (outermost)  ring 
actuator.  An  empirical  curve  fit  to  the  central  actu- 
ator profile  is  shown  as  a dotted  line  in  (a);  the  curve 
is  exp  (-3.  122  r*'  ^),  where  r is  the  radial  distance 
from  the  profile  peak  (actuator  center)  in  inches.  The 
inter-actuator  center-center  spacing  is  0.  55  in. 
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where  S is  the  interactuator  center-center  spacing,  then  Cm  = 0.  23  for  this 
RADC  beryllium  mirror. 

The  more  important  quantity,  however,  is  the  servo  coupling,  • 

5 sn 

O'Meara  has  defined  this  coupling  as  the  ratio  of  the  error  signal  induced 

in  a given  channel  by  the  displacement  of  a single  neighboring  actuator  to 

that  induced  by  a comparable  displacement  of  the  actuator  associated  with 

the  given  channel.  When  the  influence  function  is  of  the  form 


IN(r)  = exp 


(7) 


2 2 2 

where  r = x + y is  the  squared  distance  from  the  actuator  center,  then 
the  servo  coupling  is 


n2 


2/n 


'sn  2Trr(2/n)  n* 


where  n is  any  number,  T (x)  is  the  gamma  function,  and 


(8) 


I 

n 


The  quantity  p is  given  by 


p = — = (-inC  )1/n  . (10) 

r r 

o 


Reference  4 gives  the  values  of  C for  various  choices  of  C and  n.  For 
° sn  m 

the  RADC  mirror,  n = 1.  5 and  C =0.  23,  giving  C =0.  70.  For  C =0.  20, 
’ m ’ 6 6 sn  m ' 

C = 0.  67  and  if  C = 0.  15,  C =0.  58.  All  of  these  values  are  very 
sn  m sn  7 

substantial  couplings,  which  may  severely  limit  the  usable  COAT  servo 
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bandwidth  when  this  mirror  is  employed.  This  is  a point  which  we  will 
have  to  verify  experimentally. 

The  only  way  to  reduce  the  servo  coupling  once  n is  fixed,  Eq.  (3), 


is  to  narrow  the  influence  function,  i.  e.,  reduce  C . Studies  conducted 

m 


on  the  Hughes  IR&D  program,  however,  show  that  n = 2.  0 minimizes 


C for  a fixed  value  of  C . On  the  other  hand,  reducing  C will  increase 
sn  m m 


the  inter  actuator  ripple  (surface  ripple  when  all  actuators  have  equal 


excitation).  The  tradeoff  between  ripple  and  C as  a function  of  C is 

sn  m 


shown  in  Figure  13.  The  ripple  data  are  obtained  by  calculating  the  surface 
produced  by  a 99-actuator,  hexagonal-array  mirror.  Each  actuator  is 
given  a unit  ampHtude  deflection  and  all  actuators  are  assumed  to  have  the 
same  influence  function.  The  resulting  surface  is  assumed  to  be  a linear 
superposition  of  all  the  individual  actuator  influence  functions.  This  is  an 
approximation,  but  one  which  should  be  reasonably  good  as  long  as  the 
actuator  restoring  springs  dominate  the  interactuator  surface  "springs" 
caused  by  the  faceplate  stiffness  (this  is  expected  to  be  the  case  in  the  RADC 
mirror). 

Note  in  Figure  13  that  curves  (b)  and  (c)  for  the  ripple  cross  over. 
This  implies  that,  for  a given  C^,  there  is  a value  of  n in  Eq.  (7)  that 


minimizes  the  ripple.  This  optimum  exponent  is  plotted  in  Figure  14  as 


a function  of  C^.  For  large  C^,  the  value  of  n ^ approaches  2.  0,  a gaus- 


sian.  For  all  cases,  however,  nQ  ^ > 2.0,  so  the  RADC  mirror  is  expected 
to  have  a larger  surface  ripple  when  excited  than  a mirror  with  a gaussian 
influence  function. 

In  fact,  based  on  the  data  in  Figures  13  and  14,  it  appears  that  the 
influence  function  on  the  RADC  mirror  is  far  from  optimum.  A value  of  n 


between  2 and  2.  5 is  desirable,  with  Cm  ~ 10%  so  that  the  servo  coupling 


will  be  less  than  30%. 

The  experimentally  observed  surface  profile  when  two  adjacent 
actuators  are  driven  with  the  same  voltage  is  shown  in  Figure  15.  For 
reference,  the  individual  actuator  profiles  are  also  shown.  The  peak-to- 
peak  surface  ripple  is  4.  9%  of  the  surface  deflection  at  one  actuator  loca- 
tion. Figure  13  predicts  8%  ripple  if  all  actuators  are  energized.  A 
simple  superposition  of  2 actuators  with  the  same  influence  function. 
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V 


= SERVO  COUPLING  OR  PEAK-PEAK 
RIPPLE, % 


Figure  13,  Servo  coupling  C (eq.  (8))  and  peak-peak  ripple 
as  a function  of  d^Jormable  mirror  mechanical 
coupling,  Cm  (eq.  (6)).  Three  influence  function 
cases  are  shown,  a "subgaussian"  (n  = 1.  5),  a 
gaussion  (n  = 2.  0),  and  a "supergaussion"  (n  = 2.  5). 
Further  data  are  available  in  Ref.  4. 
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Figure  15.  Experimentally  observed  surface  profile  with  two 

adjacent  actuators  on  DARPA/RADC  mirror  driven 
with  equal  voltage  (=30  V rms  at  500  Hz).  The 
individual  actuator  influence  functions  that  produce 
the  resultant  surface  are  also  shown. 
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r 1.5-i 

IN  = exp  [-  (r/r  ) ’ J»  and  23%  coupling,  predicts  a 3.  3%  ripple.  The 
difference  between  this  value  and  the  observed  4.  9%  can  be  attributed  to 
the  differences  between  the  two  influence  functions  in  Figure  14;  they  are 
not  identical. 

C.  Ir radiance  Tailoring  Studies 

The  previous  contract  report*  presented  some  results  on  efforts  to 
find  a fixed  transmitter  irradiance  distribution  that  will  minimize  thermal 
blooming  effects.  The  basic  conclusion  presented  there  was  that  special  laser 
modes  and  truncated  distributions  have  little  effect  on  the  strength  of  bloom- 
ing at  a fixed  transmitted  power.  In  order  to  minimize  blooming,  the  trans- 
mitter aperture  should  be  made  as  large  as  possible  and  the  irradiance 
distribution  should  be  made  as  uniform  as  possible. 

Since  the  time  the  first  report  was  written,  we  have  found  an  error 
in  the  calculations  that  produced  Figures  7 and  10  of  that  report.  For- 
tunately, the  basic  conclusion  remains  unchanged  even  though  the  quantita- 
tive output  of  the  calculations  is  different.  The  corrected  work  has  been 
accepted  for  publication  in  the  Journal  of  the  Optical  Society  of  America, 
to  appear  in  early  1977.  Since  the  corrected  results  appear  in  that  article, 
we  present  it  here  in  its  entirety  as  an  Appendix. 
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III. 


PLANS  FOR  THE  NEXT  CONTRACT  PERIOD 


The  next  phase  of  the  contract  will  be  devoted  almost  exclusively  to 
characterizing  the  RADC  beryllium  mirror.  The  faceplate  will  be  thinned 
to  0.  125  in.  to  reduce  the  mechanical  coupling,  and  the  unit  will  be  polished 
with  19  actuators  in  place.  All  of  the  data  required  by  Amendment  No.  1 
to  the  contract  will  be  obtained.  These  data  will  form  the  substance  of  the 
3rd  interim  report  of  the  contract. 

Li  addition  to  completing  the  construction  and  characterization  of  the 

mirror,  we  will  also  complete  the  Zernike  polynomial  generator  and  set  up 

the  thermal  blooming  optical  system.  Our  present  plans  are  to  use  a 

flowing -liquid  blooming  medium  for  two  reasons:  (1)  difficulties  and  safety 

hazards  with  the  flowing-gas  cell;  (2)  ease  of  producing  substantial  blooming 

6-8 

with  the  limited  optical  power  available.  Previous  work  has  indicated  no 
signficant  differences  between  thermal  blooming  produced  in  a gaseous  or  in 
a liquid  medium. 
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PROPAGATION  OF  LASER  BEAMS  HAVING  AN  ON- AXIS  NULL 
IN  THE  PRESENCE  OF  THERMAL  BLOOMING* 

J.  E.  Pearson,  C.  Yeht,  and  W.  P.  Brown,  Jr. 

Hughes  Research  Laboratories 
3011  Malibu  Canyon  Road 
Malibu,  California  90265 

ABSTRACT 

The  propagation  of  focused  beams  having  an  on-axis  irradiance 
null  is  considered  in  the  presence  of  thermal  blooming.  Two  cases  are 
treated:  (a)  beam  profiles  that  have  an  irradiance  zero  in  the  beam  center 
at  the  focal  plane  as  well  as  the  transmitter  aperture;  (b)  beam  profiles  that 
have  an  on-axis  irradiance  null  only  at  the  transmitter.  It  is  demonstrated 
that  none  of  the  beam  profiles  considered  in  case  (a)  has  a meaningful 
advantage  over  a Gaussian  beam  profile.  Some  of  the  case  (b)  profiles  do 
produce  a larger  bloomed  irradiance  in  the  focal  plane,  particularly  when 
the  initial  intensity  distribution  is  very  uniform  in  its  non- zero  regions. 
Addition  of  a simple  central  obscuration  to  a "filled"  irradiance  distribution 
is  found  to  have  no  advantage,  however,  for  the  cases  considered. 


*To  be  published  in  the  Journal  of  the  Optical  Society  of  America.  This 
work  sponsored  in  part  by  the  Air  Force  Systems  Command's  Rome  Air 
Development  Center,  Griffiss  AFB,  NY. 

t Permanent  address:  Electrical  Engineering  Department, 

University  of  California  at  Los  Angeles 
Los  Angeles,  California  90024 
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I.  INTRODUCTION 

When  an  optical  laser  beam  is  propagated  through  an  absorbing 

medium  with  a negative  refractive  index  temperature  coefficient  (dn/dT  < 0), 

1 2 

the  phenomenon  of  thermal  blooming  ’ occurs.  In  the  presence  of  thermal 
blooming,  the  propagation  medium  acts  like  a negative  lens  that  defocuses, 
spreads,  and  distorts  the  optical  beam.  This  power -dependent  phenomenon 
can  be  particularly  harmful  when  the  optical  system  is  attempting  to  produce 
maximum  focal  plane  irradiance. 

It  has  been  suggested  that  the  use  of  a transmitted  beam  that  has  an 

annular  intensity  profile  may  provide  higher  focal  plane  irradiance  in  the 

3 

presence  of  thermal  blooming  than  that  produced  by  a Gaussian  beam.  The 
physical  reasoning  behind  this  suggestion  is  straightforward.  A Gaussian 
intensity  profile  produces  a quadratic  refractive  index  variation  in  the 
medium  that  is  minimum  at  the  beam  center  where  most  of  the  optical  pov  er 
occurs.  With  an  annular  beam  that  has  an  on-axis  null,  the  negative  lens 
induced  in  the  medium  has  a larger  refractive  index  on-axis  than  off-axis. 
The  result  is  a positive  lens  for  light  that  is  close  to  the  beam  axis,  and 
this  lens  tends  to  counteract  some  of  the  off-axis  negative  lens  beam  spread- 
ing as  the  optical  beam  propagates. 

This  article  presents  the  results  of  an  investigation  aimed  at 
quantifying  how  much  thermal  blooming  distortions  are  reduced  by  the  use 
of  laser  beam  that  have  an  on-axis  intensity  zero.  Such  beams  can  be 
produced  either  by  apodization  (truncation),  or  by  use  of  laser  resonators 
that  have  annular  gain  profiles,  ^ or  by  many  unstable  resonator  con- 
figurations such  as  those  used  in  many  high  energy  laser  devices.  The 
investigations  were  performed  by  computer  simulation  of  the  propagation 
problem.  The  computer  code^  solves  the  coupled  differential  equations 
that  describe  both  the  medium  dynamics  and  the  diffraction- propagation 
of  a focused  laser  beam.  Convection-dominated  heat  transfer  produced  by 
a transverse  wind  is  assumed  and  the  quantity  of  interest  is  the  peak 
focal  plane  irradiance. 

Two  types  of  annular  or  "on-axis-null"  intensity  profiles  are  treated. 
Section  II  discusses  those  that  have  an  intensity  zero  in  the  beam  center. 
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both  in  the  focal  plane  (far -field)  as  well  as  at  the  transmitter  (near -field) 
and  thus  are  free- space  propagation  modes.  Section  III  discusses  those 
profiles  that  have  a zero  in  the  beam  center  only  at  the  initial  transmitter 
plane.  Due  to  diffraction,  the  initial  on-axis  null  fills  in  as  these  beams 
propagate  so  that  the  on-axis  irradiance  is  not  zero  at  the  focal  plane.  An 
appropriate  choice  of  the  initial  profile  can  indeed  increase  the  maximum 
focal  plane  irradiance.  The  standard  of  comparison  in  all  cases  is  the 
focal  plane  irradiance  produced  by  a Gaussian  beam  that  is  truncated  at  the 
10%  intensity  radius. 

II.  "ON-AXIS-NULL"  IRRADIANCE  PROFILE  AT  THE  TRANSMITTING 
APERTURE  (NEAR-FIELD)  AND  AT  THE  FOCAL  PLANE 
(FAR -FIELD) 

In  this  section,  we  consider  cases  where  the  irradiance  null  in  the 

middle  of  the  beam  is  retained  as  the  focused  beam  is  propagated  to  the 

target.  To  first  demonstrate  that  an  appropriately  chosen  transmitter 

irradiance  distribution  can  indeed  be  focused  to  produce  an  on-axis -null 

focal-plane  irradiance  distribution,  we  start  with  a transmitter  aperture 
4 

field  given  by 


E(rl,ei)  = Amfm(rj)  cos(m0j),  (m  = 


0,  1,  2 


9 


(1) 


where  Am  is  the  maximum  field  magnitude,  fm(rj)  niay  be  a continuous  or 
a discontinuous  function  of  r^,  and  (r ^ ai?  the  transverse  polar  coor- 
dinates of  the  transmitting  aperture.  From  Fresnel  diffraction  formulae, 
we  obtain  the  field  u at  the  focal  plane: 


eJkz 

U (x  , y , Z)  = -rr 

O O jXz 


\ . k 2 2 1 

exp  viz  x0  + y0 1 * 
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^AO^epexp)-  jff  (xoXj  +yoYj)}  dXjdyj, 


(2) 
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where  z is  the  focal  distance,  \ the  free- space  wavelength,  k = 2n/\,  (x0»yQ) 
are  the  transverse  rectangular  coordinates  at  the  focal  plane,  and  fxj,y-j) 
are  the  transverse  rectangular  coordinates  at  the  transmitting  aperture 
plane.  Changing  Eq.  (2)  into  polar  coordinates  and  carrying  out  the 
angular  integration,  we  find 
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The  on-axis  irradiance  occurs  at  rQ  = 0 and  is  given  by 


, , f 2tt  cos(m0  ) I ‘ 

“<°.v”l2  = L — 


I2  (0)  . 

m 


(5) 


Since  J (0)  = 0 when  m 4 0,  Eq.  (4)  indicates  that  I =0  for  m 4 0.  An 
m m 

on-axis-rull  focal  plane  irradiance  is  thus  obtained  for  all  m 4 0.  When 
m = 0,  however,  the  on-axis  focal  plane  irradiance  may  take  on  any  value, 
including  zero,  depending  upon  the  initial  radial  dependence  of  u. 
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The  particular  transmitter  irradiance  profiles  that  we  have  chosen 
to  study  are  given  in  Eqs.  (6)  and  (7)  and  are  illustrated  in  Fig.  1.  The 
truncated  Gaussian  beam  we  are  using  for  comparison  purposes  is  also 
shown  in  Fig.  1 and  defined  in  Eq.  (8). 


-r2/  2 

Case  . |2  2 FP°  2 

(Il-a)  lul  = Aj  cos  6j,  for  rj  < aQ 


for  r.  > a 
1 o 


Case 


(Il-b)  |uf  - a\ 


(LliJ  _l 
V % / p2 


cos  0.,  for  b < r < a 
1 o 1 o 


= 0,  otherwise 


Truncated 

Gaussian 


2,  2 
-rl/po' 


and  b < r < a 
o m o 


for  r , < a 
1 o 


otherwise 


For  these  beams,  the  first  two  of  which  can  be  produced  by  coaxial- 

4 

resonator,  annular-gain  lasers,  pQ  is  the  e-folding  intensity  radius 
of  the  Gaussian  beam  in  Eq.  (8),  p is  a constant  less  than  or  equal  to  1, 


r is  the  adjustable  maximum  intensity  radius,  and  b and  a are  the 
inner  and  outer  radii  of  a coaxial  laser  resonator.  The  beams  in  Eqs.  (6) 
and  (7)  are  assumed  oriented  in  such  a way  that  the  6^  = 7r/2  axis  is 
parallel  to  the  transverse  wind  direction.  This  orientation  produces 
minimum  blooming  since  the  heated  air  from  one  main  transmitter  lobe 
does  not  pass  through  the  region  heated  by  the  other  main  lobe.  The  field 
amplitudes  A^,  A^,  and  A^  are  chosen  so  that  the  total  transmitted  power 
is  constant.  The  expressions  used  to  obtain  these  amplitudes  are  derived 
in  the  appendix. 

Figure  2 compares  the  initial  transmitter  irradiance  distributions 
of  the  three  beams  to  the  focal-plane  irradiance  distributions  for  a repre- 
sentative transmitter  power.  The  parameters  necessary  to  specify  the 
blooming  strength  are  given  in  the  appendix.  For  all  the  computational 
results  presented  here,  we  have  chosen  values  for  the  various  constants  to 
correspond  to  the  obscuration  ratio  in  some  laser  pointer /tracker  systems 
now  in  use,  and  to  give  the  most  uniform  initial  intensity  distribution;  we  have 
found  that  thermal  blooming  is  minimized  for  uniform  initial  intensity 
distribution;  we  have  found  that  thermal  blooming  is  minimized  for  uniform 

initial  intensity  profiles.  The  values  used  in  our  calculations  are  a /p 

o r o 

= 1.  52,  aQ/bQ  = 7,  r^/a^  = 2/3,  and  p = 1.  Note  that  in  the  absence  of 
blooming  effects,  the  focal  plane  distributions  would  look  like  the  transmitter 
distributions  for  each  case. 

The  peak  focal-plane  irradiance  as  a function  of  total  transmitted 
laser  power  for  these  three  beam  cases  is  given  in  Fig.  3.  Both  the 
irradiance  and  the  transmitter  power  are  plotted  in  normalized  units  as 
discussed  in  the  appendix.  The  maximum  focal-plane  irradiance  for 
Case  (Il-b)  is  a factor  of  1.  2 larger  than  that  obtained  with  the  truncated 
Gaussian  beam.  The  higher  irradiance  is  obtained  by  increasing  the  total 
laser  power  at  the  aperture  plane  by  a factor  of  2.  75.  Case  (Il-a)  yields 
lower  peak  focal-plane  irradiance  even  with  higher  total  laser  power. 

Hence,  no  significant  improvement  is  obtained  with  transmitting  intensity 
profiles  tailored  according  to  cases  (Il-a)  or  ( II - b ) . 
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Figure  2.  Pseudo-gray  scale  plots  of  transmitter  and  focal-plane 

irradiance  distributions  for  the  three  beams  in  Fig.  1.  A 
moderate  amount  of  thermal  blooming  is. present  in  the  pro- 
pagation path. 
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= FOCAL  PLANE  IRRADIANCE,  ARBITRARY  UNITS 


Figure  3.  Peak  focal-plane  irradiance  versus  total  transmitter  power  for 
the  same  three  cases  as  in  Figs.  1 and  2.  The  axis  units  are 
normalized  values  as  discussed  in  the  Appendix. 


HI.  "ON-AXIS-NULL"  INTENSITY  PROFILE  ONLY  AT  THE 
TRANSMITTING  APERTURE  (NEAR-FIELD) 

We  now  treat  cases  for  which  the  on-axis  intensity  null  occurs  only 
at  the  transmitting  aperture.  In  other  words,  under  thermal  blooming 
conditions,  the  cooler  region  in  the  middle  of  the  beam  does  not  exist 
throughout  the  focused  propagation  path  to  the  focal  plane.  Four  specific 
aperture  intensity  profiles  are  considered  and  compared  to  the  truncated 
Gaussian  in  Eq.  (8): 
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The  beam  amplitudes  A3,  A 4,  Ag,  A^  are  chosen  to  keep  the  total 
transmitter  power  constant.  The  calculations  for  these  amplitudes  are 
given  in  the  appendix.  Computer-generated  plots  of  the  initial  beam 
pforiles  are  illustrated  in  Fig.  4. 

The  transmitter  aperture  intensity  profiles  as  well  as  the  corre- 
sponding focal  plane  intensity  profiles  for  these  cases  are  displayed  in 
Fig.  5 for  a representative  laser  power  level.  Figure  6 compares  the 
normalized  peak  focal  pline  irradiance  as  a function  of  normalized  total 
input  transmitting  power  for  these  four  cases  to  that  for  the  truncated 
Gaussian  case.  Unlike  the  cases  treated  in  Section  II,  significantly  larger 
peak  focal  plane  irradiance  is  achieved  for  two  of  the  four  cases  considered 
as  compared  with  that  for  the  truncated  solid  Gaussian.  Case  (Ill-d)  yields 
a peak  focal  plane  irradiance  which  is  the  same  as  that  for  the  truncated 
Gaussian  beam. 

The  intensity  profile  for  Case  (Ill-a)  is  obtained  by  combining  the 

TEMnn  and  TEM*  modes  of  a cylindrical  cavity  oscillating  in  phase 
U U 1 u ^ 

opposition  (Laguerre -Gaussian  modes  ).  The  peak  focal  plane  irradiance 
achieved  with  this  profile  is  almost  twice  that  of  the  truncated  solid  Gaussian 
beam  with  a total  laser  power  roughly  2.  5 times  that  for  the  truncated  solid 
Gaussian  beam.  Other  ways  to  increase  focal  plane  irradiance  by  this 

g 

amount  such  as  coherent  optical  adaptive  techniques  (COAT)  may  be  more 
practical  and  efficient  if  this  type  of  mode  is  difficult  to  produce  in  practice 
in  high  power  lasers. 

In  cases  (IU-b)  and  (III-c),  the  peak  focal  plane  irradiance  is  also 
about  twice  the  value  for  the  truncated  solid  Gaussian  beam,  but  as  shown 
in  Fig.  6,  this  increase  is  obtained  at  about  90%  of  the  total  input  laser 
power  of  the  truncated  solid  Gaussian  beam.  A very  uniform  intensity 
distribution  is  thus  desirable  for  reducing  blooming,  but  may  be  very 
difficult  to  obtain  in  practice. 

IV.  CONCLUSION 

These  results  lead  us  to  conclude  that  many  special  laser  modes  or 
intensity  distributions  that  have  an  on-axis  null  or  minimum  do  not  offer 
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4629  5 


CASE  nr -a 


case  m-b 


Figure  4.  Three-dimensional  plots  of  transmitter  irradiance  profiles  that 
have  an  on-axis  null  initially  that  fills  in  as  the  beams  pro- 
pagate. Case  Ill-a  is  a mixture  of  cylindrical  modes  TEMqq 
and  TEMqj  oscillating  in  phase  opposition.  Case  Ill-b  is 
similar  to  Case  II-bf  except  | u | 2 has  no  0 -dependence. 

Case  III-c  is  a uniform  annular  beam.  Case  Ill-d  is  an 
annularly-truncated  Gaussian  beam. 
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Figure  5.  Pseudo-gray  scale  plots  of  transmitter  and  focal-plane 

irradiance  distributions  for  the  three  beams  in  Fig.  4.  A 
moderate  amount  of  thermal  blooming  is  present  in  the 
propagation  path. 
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FOCAL  PLANE  IRRADIANCE,  ARBITRARY  UNITS 


Figure  6.  Peak  focal-plane  irradiance  versus  total  transmitter  power 
for  the  four  beams  in  Figs.  5 and  6.  The  axis  units  are 
normalized  as  discussed  in  the  Appendix. 
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any  significant  advantages  for  reducing  thermal  blooming  when  the  goal  is 
to  achieve  maximum  focal-plane  irradiance  for  a given  transmitter  power. 

In  particular,  the  laser  modes  considered  here  that  produce  an  on-axis  null 
in  both  the  near-field  and  far-field  have  no  advantages  over  intensity  dis- 
tributions that  have  the  null  only  in  the  near -field  of  the  transmitter. 
Significantly  less  blooming  does  occur,  however,  for  certain  annular,  cir- 
cularly symmetric  laser  modes  and  for  initial  irradiance  profiles  that  are 
very  uniform  in  the  non- zero  regions  (a  uniform  annulus,  for  example). 

One  question  that  is  frequently  asked  is  whether  the  "holes"  intro- 
duced into  beams  by  unstable  resonators  should  be  eliminated  or  whether 
they  may  provide  some  advantage.  Although  we  have  not  considered  a 
sufficient  number  of  cases  to  make  a definite  conclusion,  our  studies  to 
date  of  filled  beams  and  of  beams  with  a central  hole  or  obscuration  indicate 
that  the  hole  makes  very  little  difference,  at  least  as  far  as  thermal  bloom- 
ing is  concerned.  When  the  goal  is  to  minimize  thermal  blooming  dis- 
tortions so  that  maximum  focal-plane  irradiance  is  achieved,  then  it  is  much 
more  profitable  to  force  the  initial  beam  intensity  profile  to  be  as  uniform  as 
possible  in  its  non- zero  regions  rather  than  attempt  to  produce  or  eliminate 
central  holes  in  the  beam.  Increasing  the  initial  transmitted  beam  diameter 

is  also  quite  effective  in  reducing  blooming  distortions,  particularly  when 

o 

coupled  with  adaptive  optical  techniques. 

Since  it  may  be  impractical  in  many  high  power  laser  applications  to 
achieve  either  sufficiently  large  apertures  or  to  produce  the  uniform  or 

special  intensity  distributions  that  minimize  thermal  blooming,  other 

1011  812 
irradiance-tailoring  techniques  ' or  phasefront-tailoring  techniques  ’ 

may  be  more  useful.  Additional  work  is  now  in  progress  in  both  of  these 

areas. 
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APPENDIX 


In  this  appendix,  we  derive  the  normalizations  used  in  the  text  for 
the  focal-plane  irradiance  and  for  the  field  amplitudes  for  the  various  beams 
in  terms  of  the  total  laser  transmitter  power.  The  total  transmitted  power 
is  given  by 


lE 1 2 dAT  ■ ^rrr  4 AmFi«V%>  <A-‘> 

where  A^,  is  the  aperture  area,  E is  the  field  amplitude  a s defined  in 
Eq.  (1),  and  is  a constant  given  by 

pa  /\TZ  p 
j o ro 

Fl<a0/p0>  = I fm(N^2  PoX)  xdx  . (A-2) 

^ 0 

The  maximum  aperture  dimension  is  aQ  and  pQ  is  the  e-folding  intensity 
radius  of  the  Gaussian  beam  in  Eq.  (8). 

The  focal  plane  irradiance  as  plotted  in  Figs.  3 and  6 is  found 
from  Eqs.  (3)  and  (4): 
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with  Im  defined  by  Eq.  (4).  We  are  interested  in  the  maximum  focal-plane 

irradiance,  which  will  occur  at  0 =0  and  at  some  r , depending  on  the 

o o 

initial  irradiance  distribution.  Using  0q  = 0 and  combining  Eqs.  (A-2), 
(A-3),  and  Eq.  (4),  the  irradiance  can  be  written  as 
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u ^ = 4t r(m  + 1 ) 
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(A-4) 


where  a is  the  propagation  medium  absorption,  Fj  is  defined  in  Eq.  (A-2), 
and  is  given  by 
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and 
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We  have  evaluated  the  integral  function  F^  in  Eq.  (A-5)  for  each  of 

the  7 radial  f functions  discussed  in  the  text  (Eqs.  (6)-(12))  to  find  the 
m 

maximum  value  of  F^  as  a function  of  the  normalized  parameter  p.  The 

values  a /p  =1.  52,  a /b  - 7,  r /a  = 2/3,  and  p = 1 were  used  for  this 
o o o o m o 

evaluation.  Table  A-l  gives  the  results  of  this  analysis.  To  obtain 
numerical  results,  we  chose  \ = 10.  6 pm,  z = 2 km,  and  pQ  = 23  cm 
(corresponding  to  a = 35  cm).  The  irradiance  plotted  as  the  ordinate 
axes  in  Figs.  (3)  and  (6)  is  thus  given  by 
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in  units  of  power  per  square  centimeter,  with  (S.R.  ) being  the  strehl 
ratio  determined  by  the  thermal  blooming  computer  code  (S.R.  = 1 for 
free- space  propagation).  For  the  blooming  calculations,  the  following 
parameters  were  taken: 
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absorption  =2x10  m 
transverse  wind  velocity  =10  m/sec 
beam  slew  rate  = 20  mrad/sec 


Table  A-l.  Values  for  Peak  Focal- Plane  Irradiance 
Calculations 


Case 

m 

F1 

(F^)  max. 

(p)  max. 

Inf.  Gau  s . 

0 

0.  25 

0.  50 

0 

Trunc.  Gaus. 

0 

0.  225 

0.  342 

0 

Il-a 

1 

0.  225 

0.  167 

2.  50 

H-b 

l 

0.  508 

0.  279 

2.  31 

m-a 

0 

0. 0509 

0.  1 6 1 

0 

Ill-b 

0 

0.  508 

0.  536 

0 

HI-c 

0 

0.  566 

0.  566 

0 

Hl-d 

0 

0.  214 

0.  331 

0 
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